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1     INTRODUCTION 


Abstract 
We  present  a  self-contained  expKieition  of  the  theory  of  Grobner 
basis  and  its  applications. 


1      Introduction 

These  notes  attempt  to  present  in  a  self-contained  manner  the  basic  facts 
about  the  theory  of  Grobner  basis  and  related  algorithms.  Except  for  an 
excellent  survey  on  the  subject  by  Buchberger  [Buchberger  1985]  the  liter- 
ature on  the  subject  is  somewhat  scattered.  Since  our  interest  is  in  com- 
plexity of  algorithms,  we  attempt  to  give  quantitative  forms  of  theorems 
whenever  possible.  Of  course,  much  of  this  material  is  ultimately  owed  to 
B.  Buchberger  even  if  not  explicitly  mentioned. 

A  Grobner  basis  is  a  special  basis  for  a  multivariate  polynomial  ideal 
over  a  field  with  certain  attractive  computational  properties.  It  turns  out 
that  many  important  computational  problems  involving  ideals  can  be  eaisily 
solved  once  we  have  a  Grobner  basis  for  the  ideals.  Such  bases  were  first 
defined  by  Hironaka  in  1964  who  called  them  standard  bases.  Buchberger 
independently  defined  the  concept  in  his  PhD  thesis  in  1965,  naming  it  in 
honor  of  his  teacher,  W.  Grobner.  Hironaka  only  proved  the  existence  of 
such  bases;  Buchberger  gave  an  aJgorithm  (the  'basic  algorithm'  in  section 
3)  to  construct  them. 

More  generally,  this  area  falls  under  that  of  computational  algebraic  ge- 
ometry. Hilbert  is  often  cited  as  launching  the  niainstream  modern  math- 
ematical trend  which  favors  non-constructive  results.  In  particular  Hilbert 
solved  the  outstanding  problem  in  the  theory  of  invariants  via  his  cele- 
brated Basis  Theorem:  that  every  polynomial  ideal  has  a  finite  bjisis.  In  a 
sense,  he  invented  his  solution  of  Basis  Theorem  because  his  solution  was 
not  in  the  spirit  of  what  the  geometers  of  his  day  (Gordon  wzis  the  princi- 
pal representative)  were  trying  to  do.  The  geometers  were  trying  to  give 
constructive  or  algorithmic  solutions.  Hilbert  tiirned  it  into  an  existence 
question,  although  he  Weis  later  able  to  return  to  the  original  problem  to 
give  a  satisfactory  construction.  As  for  the  trend  towards  non-constructive 
abstractions,  the  subject  of  algebraic  geometry  itself  is  the  best  illustration 
of  the  phenomenon  (see  the  historical  survey  of  the  whole  subject  till  present 
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day  in  [Dieudonn^  1985]).  Meanwhile,  the  constructive  spirit  in  algebraic 
geometry  has  survived,  though  not  thrived.  A  key  paper  in  this  regards  is 
that  of  Hermann  [Hermann  1926].  Another  more  modem  attempt  to  re- 
turn to  the  classical  questions  is  Seidenberg  [Seidenberg  1974].  Since  then  a 
number  of  papers  have  begun  to  revitalize  this  area,  and  among  the  reasons 
for  this  is  the  promise  of  practical  computational  algebra  systems. 

We  would  like  to  noake  two  remarks  here:  a  prominent  algebraic  geome- 
ter advocating  the  constructive  viewpoint  is  Abhyankar  (see  his  delightful 
historical  expos6  in  [Abhyankar  1976]).  As  computer  scientists,  we  take 
Professor  Abhyankar's  viewpoint  to  the  extreme':  we  regard  the  existence 
of  a  construction  only  as  a  first  step  towards  a  precise  clzissification  of  the 
inherent  computational  complexity  of  a  computational  algebraic  problem. 
Studies  motivated  by  this  extreme  requirement  might  be  termed  algorith- 
mic algebraic  geometry.  It  is  our  belief  that  this  standpoint  can  transform 
a  large  area  of  classical  algebraic  geometry  to  a  new  level  of  beauty  and 
understanding.  Already  this  algorithmic  standpoint  has  revived  the  study 
of  elementary  Euclidean  geometry.  The  other  remark  is  this:  the  study  of 
algorithmic  algebraic  geometry  Implies  quantitative  algebraic  geometry  in 
which  a  variety  of  bounds  are  sought.  To  illustrate  this  p>oint,  Hilbert's 
baisis  theorem  tells  us  that  all  ascending  chains  of  polynomial  ideals  are 
finite.  We  seek  to  bound  the  length  of  such  chains.  (It  appears  that  we 
cannot  bound  the  length  as  a  fimction  of  the  initial  ideal  alone).  Or  again, 
Hillbert's  Nullstellensatz  says  that  if  /  vanishes  at  all  the  zeroes  of  am  ideal 
/  C  ik[ii, . . . ,  i„]  where  k  is  algebraically  closed  then  f^  £  I  for  some  d.  In 
fact  it  had  been  noted  that  d  depends  only  on  /,  not  on  /.  We  again  ask 
for  good  bounds  on  d  as  a  function  of  the  size  of  /. 


2      Basic  Terminology 

Let  us  fix  a  polynomial  ring  R  =  K[xi, . . . ,  x„],  n  >  1,  and  K  is  any  field. 
Given  elements  a,b  e  Rv/e  say  that  a  divides  b,  equivalently  fc  is  a  multiple 
of  a,  if  6  =  ac  for  some  c  E  R.  If  Oi, . . . ,  a„  are  elements  in  R  and  F  C  R 


^Alternatively,  as  computer  scientists,  we  have  more  precise  notions  of  constructiveness 
and  a  rich  vocabulary  to  make  finer  distinctions. 
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then  the  ideals  generated  by  {ai, . . .  ,o„}  and  F  (respectively)  are  denoted 
by  (ai,...,a„)  and  {F). 

A  power  product  is  an  element  of  R  of  the  form 

p  =  i^'i'j'-i'„",  c.>0. 

Sometimes  power  products  are  also  called  terms.  The  total  degree  of  p  is 
deg[p)  =  52"^  1  c,.  The  degree  of  p  in  any  variable  i,  is  (feff,_(p)  =  c^.  The 
/eas<  common  multiple  (LCM)  of  two  power  products  p  =  xj'xj'  " " "  ^n"  and 
9  =  xf'i^---i^  isgivenbyxr*^'''"^---ir^t'^''"^  A  monomto/ is  a  term 
of  the  form  ap  where  a  ^  K  and  p  is  a  power  product.  The  length  of  a 
polynomial  is  the  number  of  monomials  that  occurs  in  it.  The  LCM  of  two 
monomials  am  and  a'm',  where  a,a'  e  K  and  m,m'  are  power  products,  is 
defined  as  follows: 

LCM{am,am')  =  aa'  •  LCM{m,m'). 

We  remark  that  if  K  is  an  Euclidean  domain  then  it  is  better  to  define  the 
LCM  as  LCM{a,a')  ■  LCM{m,m'). 

Let  PP  =  PP(ii,. ..  ,i„)  be  the  set  of  all  power  products  involving 
Zi, . . . ,  z„.   A  total  ordering  <  on  the  set  of  power  products  is  said  to  be 

A 

admissible  if  for  all  power  products  p,  q,  (i)  1  <  p  for  all  p  €  PP  and  (ii) 

(Va  e  PP)  p  <  9  =^  ap<aq 

The  usual  examples  of  admissible  orderings  are  lexicographical  ordering  etnd 
the  total  degree  ordering  in  which  p<  9  if  either  deg{p)  <  deg{g)  or  else,  in 

A 

case  deg{p)  =  deg{g),  then  p  is  lexicographically  less  than  g.  Note  that 
both  the  lexicographical  aind  total  degree  ordering  are  completely  specified 
by  a  total  ordering  on  the  variables,  say,  Xi  <  •  •  •  <  x„. 

A  A 

Remark:  Admissible  orderings  are  sometimes  called  term  orderings  in 
the  literature.  Any  ordering  that  satisfies  (ii)  (but  not  necessarily  (i)) 
may  be  called  semi-admissible.  For  instance,  the  reverse  of  any  admissible 
ordering  is  semi-admissible.  Lazard  has  pointed  out  that  the  reverse  lex- 
icographical ordering  has  interesting  algorithmic  properties  [Lazard].   One 


2    BASIC  TERMINOLOGY 


crude  classification  of  all  semi-admissible  orderings  is  given  by  the  total 
order  imposed  on  the  set  {l,  ij,  ij, . . . ,  x„} 

Henceforth,  fix  <  to  be  any  admissible  ordering.  We  also  write  f  >  g 
if  f  ^  g  and  g<  f.    Let  /  e  J?  be  a  polynomial.    The  head  monomial 

A 

Hmono(/)  of  /  is  the  monomial  in  /  whose  power  product  is  largest  relative 
to  <.  (If  /  =  0  we  define  Hmono(/)  =  0.)  For  instance,  relative  to  the  total 

degree  ordering,  the  head  term  of  f  —  Axy  +  y  —  5  is  n{f)  =  4xy.  We  may 
define  the  head  term  Hterm(/)  to  be  the  power  product  in  Hmono(/);  the 
definition  of  head  coefficient  Hcoef(/)  of  /  is  then  clear.  Thus  Hmono(/)  = 
Hcoef(/)  •  Hterm(/).  Ocassionally,  it  is  useful  to  have  the  notation  Tail(/) 
for  /  —  Hmono(/). 

A  basic  concept  in  Grobner  basis  is  the  idea  of  reduction.  Given  two 
polynomials  f,g  £  R,  we  say  /  is  reducible  by  g  if  Hmono(ff)  divides  some 
monomial  m  in  /.  Say  m  =  c  •  Hmono(y).  Then  we  say  the  polynomial 
h  =  f  —  c-  g  is  the  reduct  of  /  by  </  and  denote  the  relationship  by 

f-^h. 

We  say  that  the  monomial  m  (or  the  corresponding  power  product  p)  is 
eliminated  by  application  of  g  in  this  case.  If  G  is  a  set  of  polynomials,  we 
write  /  — ►  h  if  f  — ^  h  holds  for  some  g  E  G.  U  there  is  a  finite  sequence 
hi,h2,. .  .,hn  (n  >  1)  such  that  hi  =  f,hn  =  h  and  hi — ^^+1  for  t  = 
1,. . .  ,n  —  1,  then  we  write 

f^h. 

If  /  is  not  reducible  by  any  </  G  G,  we  indicate  this  by  writing  /  — >  /.  We 
say  /i  is  a  G-normal  form  of  f  if  f  — >  h  — >  h,  and  we  write  NFg(/)  for  the 
set  of  all  G-normal  forms  of  /.  It  is  important  to  realize  that  the  G-normal 
form  of  /  is  not  unique  in  general,  and  the  central  idea  in  Grobner  basis  is 
to  enlarge  G  so  that  it  becomes  unique.  Finally,  we  are  ready  for  the  main 
definition:  A  finite  set  G  C  i?  is  seiid  to  be  a  Grobner  basis  (for  the  ideal 
generated  by  G)  if  the  G-normal  form  of  every  polynomial  /  is  unique,  i.e., 
|NFa(/)|  =  l. 

Remark:  Computationally,  it  may  be  more  efficient  (asymptotically) 
to  use  the  reflexive  transitive  closure  of  — >.    If  =$•  denotes  this  closure, 


3    NORMAL  FORM  ALGORITHM  6 

then  f  ^  h  implies  that  no  monomial  in  /i  is  a  multiple  of  Hmono(3). 

Exercise.  Prove  that  h  is  imiquely  determined  by  /  and  g. 

It  turns  out  that  the  concept  of  a  Grobner  basis  is  intimately  related 

to  that  of  the  S- polynomial  S{f,g)  of  two  polynomials  f,g  £  R.    This  is 

defined  by 

m  "^ 


Hmono(/)         Hmono(^) 

where  m  =  LCM(Hmono(/),Hmono(y)).  For  example,  assuming  a  total 
degree  ordering, 

5(2i'y  +  zy  +  y'  -  3,5y'  -  y  +  i)  =  5y(iy  +  y^  -  3)  -  2i'(-y  +  i). 

Example:  Admissible  orderings  can  be  quite  complex.  The  following 
example  is  due  to  T.  Dube:  it  shows  an  ordering  on  PP(x,y,2)  with  the 
property  that  the  restriction  of  the  ordering  to  PP(x,y)  is  a  total  degree 
ordering  with  y  >  x,  x>  z*  and  y  >  z*  for  all  i  and  yet  xz  >  y.  The  rule  is 

A  A  A  A 

this.  Suppose  m  =  x'y^'z'  and  m'  =  z"  y*'  and  z**.  Then  m>  m'  iff  one  of 

the  following  sequence  of  tests  yields  an  affirmative  answer  when  applied 
in  the  indicated  order:  (1)  a  +  b>  a'  +  b',  (2)  c  >  c\  (3)  a  <  a!.  One  should 
check  that  this  is  an  admissible  ordering  and  has  the  claimed  properties. 


3     Normal  Form  Algorithm 

Henceforth  we  fix  some  admissible  ordering  <  on  monomials.  For  any  finite 

A 

set  F  C  R,  we  define  a  "non-deterministic"  aJgorithm  which  for  any  input 
polynomial  /,  computes  an  element  in  NFf  (/).  The  algorithm  is  trivial: 
apply  any  reduction  — >,  y  E  F,  to  transform  the  input  polynomizil  /.  Now 
as  long  as  the  (transformed  input)  polynomial  is  not  in  NF/'(/),  repeat  the 
reduction.  We  write  nfj?(/)  for  the  polynomial  produced  by  this  aJgorithm. 
All  reductions  in  this  section  are  relative  to  a  fixed  but  arbitrary  set  F. 

Before  proving  the  termination  of  this  algorithm  we  prove  Dixon's  lemma 
[Dixon  1913]  and  two  general  results  on  well-ordering. 

Lemma  1  [Dixon]  Every  stt  X  C  PP  of  monomials  contains  a  finite  subset 
F  C  X  such  that  each  m  ^  X  is  a  multiple  of  some  monomial  in  F. 
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Proof.  We  use  induction  on  the  number  n  of  variables.  If  n  =  1  then  we  let 
F  consist  of  the  unique  polynomial  in  X  of  minimum  degree.  So  we  may 
assume  n  >  1.  Pick  any  /o  €  X  and  say 


Then  every  m  G  A"  that  is  not  divisible  by  /o  belongs  to  one  of  J2i=i  c» 
different  sets:  let  t  =  l,...,n  and  v  =  0, 1, . . .  ,Cj  —  1.  Then  the  set  X,_„ 
consists  of  those  monomials  m  G  X  such  that  degi_(m)  =  v.  Let  X,'„  denote 
the  set  of  monomials  obtained  by  omitting  the  factor  x^  from  monomials 
in  X,  „.  By  inductive  hypothesis,  there  exists  finite  subsets  /"/^  C  X,'„  such 
that  each  monomiaJ  in  X,'„  is  a  multiple  of  some  monomial  in  F/^.  We 
obtain  Fj  „  as  {m  ■  x"  :  m  G  F/^}.  It  is  then  clear  that  every  monomial  in 
X  is  a  multiple  of  some  monomial  in  the  finite  set 

{/o}uUn- 
Q.E.D. 
Lemma  2  Every  admissible  ordering  <  on  PP  is  a  well-ordering. 

A 

Proof.  This  is  an  easy  consequence  of  the  Dixon's  lemma.  Suppose  we  have 
an  infinite  descending  sequence  of  monomials 

mi  >  mj  >  •  •  •  >  m,  >  •  •  • . 

A  A  A  A 

Let  X  =  {'Tii,»7i2, . . .  ,m,, . . .}  and  let  F  C  X  be  a  finite  subset  such  that 
every  m  G  X  is  a  multiple  of  some  monomial  in  F.  Let  m'  be  the  mono- 
mial that  is  smallest  in  F  imder  the  ordering  <.    Since  every  m  G  X  is 

A 

greater  than  some  monomial  in  F,  we  have  m'  <  m.   Thus  m'  terminates 

A 

the  descending  sequence,  contradicting  our  assumption  that  the  sequence 
is  infinite.  Q.E.D. 

For  the  next  result,  we  need  a  definition.  Let  X  be  any  set  with  a  total 
ordering  <'  and  let  S[X)  be  the  set  of  all  finite  decreasing  sequences  of 
elements  of  X: 

5(X)  =  {(ii,X2,...,i„)    :   X,  G  X,ii  >  X2  >  •••  >  !„} 
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Let  S{X)  have  the  following  induced  total-ordering: 

(ii,z,,...,x„)  <'  (yi,y2,...,ym) 

if  either  for  some  t  <  rain(n, m),  Xi  =  yi,...,z,  =  y,  and  i^+i  <  y,+  i,  or 
else  the  sequence  (ij , . . . ,  i„)  is  a  prefix  of  the  sequence  (yi, . . . ,  ym)  (thus 
n  <  m). 

Lemma  Z  If  X  is  well-ordered  by  <'  then  S{X)  is  well-ordered  under  the 
induced  ordering. 

Proof.  For  the  sake  of  contradiction,  suppose  Oi  >'  02  >'  •  ■  •  is  an  infinite 
descending  chain  in  S{X).  Lee  a,  =  (2^1,1,2, ...  ,2:, ,„(,)).  There  are  two 
cases. 

(i)  The  n(t)'s  are  bounded,  say  k  =  majc{n(t)  :  t  =  1,2,---}.  We  use 
induction  on  k.  We  get  an  immediate  contradiction  for  k  =  1,  so  assume 
k  >  1.  li  there  are  infinitely  many  t's  such  that  n(i)  =  1  then  we  get  a 
contradiction  from  the  subsequence  consisting  of  such  a,'s.  Hence  we  may 
assume  that  the  n(t)'8  are  all  greater  than  1.  Now  there  is  an  to  such 
that  for  all  t  >  to»  ^i,i  =  ^i+i,i-  Let  a,'  =  (2:i,j,i,,s5  •  •  •  >a^i,n(i))  he  obtained 
from  <7,  by  omitting  the  leading  item  in  the  sequence.  Then  the  sequence 
'^o'^'^o+i' *  ■  ■  constitute  a  strictly  decreasing  decreasing  infinite  chain  with 
each  <7^  of  length  <  k.  This  contradicts  the  inductive  hypothesis. 

(ii)  The  n(i)'s  are  unbounded.  By  taking  a  subsequence  if  necessary, 
we  may  assume  that  n{»)  is  strictly  increasing  in  t.  Define  m{l)  to  be  the 
largest  index  such  that  im(i).i  =  Xji  for  all  j  >  m(l).  For  each  t  >  1  define 
m(i)  to  be  the  largest  index  greater  than  Tn[i  —  1)  such  that  1,0(1), 1  =  ^y.t 
for  all  j  >  m{i).  Note  that  the  sequence 

Xm(l).\,  Xm(7).7,  Xn(S),i,  *  *  * 

is  strictly  decreasing.  This  contradicts  the  well-foundedness  of  X.  Q.E.D. 


Theorem  4   The  normal-form  algorithm  terminates. 

Proof.  We  map  a  polynomial  g  to  the  sequence  of  monomials  g  =  (mi, . . . ,  m^t) 
where  m,  are  the  monomials  occurring  in  g  and  mi  >  mj  >  •  •  •  >  mt.  By  the 
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last  two  lemmas,  the  set  of  p's  are  well-ordered  under  the  induced  ordering 
<•  It  is  seen  that  if  g — »/i  then  ^>h.  The  termination  of  the  algorithm 

A  ^ 

is  equivalent  to  the  well-fovmdedness  of  the  induced  ordering.  Q.E.D. 


4      Bounds  on  Normal  Form  Algorithms 

We  now  quantify  the  termination  process  of  the  normal  algorithm.  More 
precisely,  let  F  C  i2  be  a  finite  set  and  g  E  R.  Let 

g  =  go — ^ffi — ' ^fft  =  nff(?)-  (1) 

Our  goal  is  to  bound  the  maximum  value  of  k  in  the  reduction  sequence 
(1).  It  is  not  surprising  that  the  bounds  will  depend  on  (a)  the  strategy  for 
choosing  reduction  steps  and  (b)  the  choice  of  admissible  ordering.  Note 
that  the  strategy  in  (a)  includes  the  one  of  'no  strategy'  (as  implied  in 
the  description  of  the  normal  form  algorithm  in  the  last  section) .  We  will 
consider  two  important  admissible  orderings:  total  degree  ordering  and 
lexicographical  ordering.  The  following  notations  will  be  used 

F  =  a.  finite  subset  of  R 

g  =  input  polynomial 

n  =  the  number  of  variables,  R  =  K[xi,. . .  ,Xn] 

d  =  the  maximum  degree  in  any  variable  x,  of  polynomials  in  F 

t  =  the  maximum  length  of  any  polynomial  in  F 

D  =  the  maximum  degree  in  any  variable  x,  of  g 

m  =  the  number  of  polynomials  in  F. 

First  we  consider  the  easier  case  of  total  degree  ordering. 

Lemma  5  //  the  total  degree  ordering  is  used,  the  number  of  reductions 
steps  to  derive  nf/?(g)  starting  from  g  is  bounded  by  2^'^^^^" . 

Proof.  Consider  the  polynomials  (/,,  (i  =  0, 1,---)  in  the  above  reduction 
sequence  (l).  The  number  of  power  products  where  the  maximima  degree 
in  any  variable  is  D  is  (D  +  1)".  Thus  the  nimaber  of  distinct  polynomials 
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with  such  power  products  (modulo  the  map  g  >-*  g)  is  2^^"*"^^".  Since  all  the 
gi  are  distinct,  the  desired  boimd  in  the  lemma  follows.  Q.E.D. 

We  next  give  a  better  bound  under  the  assumption  that  the  normal 
form  algorithm  always  chooses  to  eliminate  the  <-largest  monomial  that 

A 

could  be  eliminated.  More  precisely,  suppose  that  in  the  reduction  step 
g,-i  — »^  <7,  in  (l)  the  monomial  m,  is  eliminated.  We  say  that  the  reduction 
sequence  (1)  is  ordered  if 

mi  >  mj  >  •  •  •>  rrik. 

AAA 

In  the  sequence  (l),  assume  m,  is  eliminated  by  application  of  /,  G  F.  Thus 

9i  -  9i-i  -  OLtfi 

for  some  monomial  a,. 

Lemma  6  Under  the  total  ordering,  the  length  of  an  ordered  reduction  se- 
quence is  at  most  {D  +  l)". 

The  proof  is  immediate.  However  we  note  that  there  may  be  members 
of  NFj?(<7)  not  reachable  by  ordered  reduction  sequences  and  furthermore, 
the  shortest  reduction  sequence  from  g  to  some  /  by  an  ordered  reduc- 
tion sequence  may  actually  be  longer  than  allowing  arbitrary  reduction 
sequences. 

Now  we  consider  the  caise  where  <  is  the  lexicographical  ordering.  We 

A 

had  noted  before  that  every  admissible  ordering  can  be  crudely  classified 
by  considering  the  order  on  the  n  variables.  We  will  assume  in  the  rest  of 
this  paper  that 

a:r,>a;n-i  >  •••>Xi. 

A  A  A 

We  define  a  weighting  function  Wp  :  PP  — >  N  (N  is  the  set  of  natural 
niimbers)  as  follows: 

V^'^(xl' •••i;-)  =  ci((f  +  1)°  + ej(d+ 1)^  +  •■•  +  e„((f  +  1)"-' 

If  /  is  any  polynomial,  we  extend  the  weight  function  to  let  Wf{f)  denote 
the  weight  function  applied  to  the  heaA  term  in  /.  Note:  Unlike  the  case 
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of  total  ordering,  oiir  bounds  for  lexicographical  ordering  will  depend  on 
the  set  F.  This  is  reflected  in  the  appearance  of  the  subscript  F  in  the 
weighting  function  Wp  (and  also  Wp  to  be  introduced  below). 

Lemma  7  Let  f  ^  F  and  let  m  be  any  monomial.  Write  f  as  a  sum  of 
monomials,  f  =  /i  +  /z  +  •  •  •  +  A  where  /i  >  /:  >  •  •  •  >  A-  Then  for  any 
monomial  m,  we  have  that  WF{mfj-i)  >  Wp{mfj)  for  j  =  2, . . . ,  fc.. 

Proof.  Let  /y_i  =  xf '  •  •  •  ij|"  and  /,  =  xV  ■  ■  ■  x'n-  Then  /y_i  >  fj  iff  there 
is  some  t  =  l,...,n  such  that  d„  =  c„,(f„_i  =  e„_i, . . .  ,rf,+  i  =  e.+i  and 
di  >  e,.  We  get 

l^f(m/;_i)-VVf(m/y)     =    Wp{f,_,)-Wp{f,) 

(c,_,  -  <i._j)(d  +  1)'-'+  •••  +  (ei  -  d,){d  +  in 

We  note  that  e^  —  dj  is  at  most  d  (the  maximum  degree  of  any  polynomial 
in  F;  hence  the  summation  in  the  square  brackets  above  is  at  most 

d{d  +  1)-'  +  •  •  •  +  d(d  +  1)°  =  (<f  +  1)'-'  -  1. 

Thus  Wp{mfj.i)  >  Wpimfj).  Q.E.D. 

Now  we  define  a  weight  function  on  polynomials.  Recall  that  the  length 
of  each  polynomial  in  F  is  at  most  t.  For  any  polynomial  h  where  h  = 
hi  -\-  h2  +  ■  ■  ■  +  ht,  hi  >  hi  >  •  ■  •  >  hk,  v/e  define  its  weight  to  be 

»=i 

Note  that  if  the  length  of  /r  is  L  and  the  mzLximimi  degree  of  any  variable 
is  D  then  Wp{h,)  <  {d+l)"{D/d]  and 

Lemma  8  If  f -^  g  then  Wpif)  >  Wpig). 
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Proof.  Suppose  g  =  f  —  rnh  for  some  h  E  F  and  monomial  m.  Let  the 
monomials  of  h  be  ordered  as  /ii  >  hj  >  •  •  •.  Consider  how  the  weight  of  g 
differ  from  that  of  /.  The  decrease  in  the  weight  of  g  over  /  due  to  the 
elimination  of  m/ij  is  at  least  ^'*'''("^').  Clearly,  all  the  new  (i.e.  not  in  /) 
monomials  in  g  is  of  the  form  m/i,  for  some  t  >  1.  Since  there  are  at  most 
£  —  1  new  monomials  in  g,  the  increase  in  the  weight  of  g  over  /  due  to 
their  presence  is  at  most  (£  -  i)t^''i"^^).  Hence 

Wjr{f)  -  WF[g)     >    l^"^"^^)  -{t-  i)£W('^^) 

_      fWr{mhi)  _  fl  +  Wr(mh3)  _^_  fW(mh2) 

By  the  previous  lemma,  Wplmhi)  >  Wf  (m/ij).  Hence  Wplf]  -  Wpig)  > 
fWimh,)  >  0.  Q.E.D. 

We  immediately  conclude: 

Theorem  9  Assuming  the  lexicographical  ordering,  the  length  of  any  se- 
quence of  reductions  beginning  from  an  input  polynomial  g  is  at  most 

where  L  is  the  length  of  g,  I  is  the  maximum  length  of  a  polynomial  in  F 
and  D  is  the  maximum  degree  of  g  in  any  variable. 

In  a  subsequent  paper  [Dub^,  Mishra  and  Yap  1986],  we  will  extend  the 
preceding  analysis  of  normaJ  form  algorithms  to  arbitrary  admissible  order- 
ings. 

5      Characterizations  of  Grobner  Basis 

Recall  that  a  finite  set  G  C  fi  is  a  Grobner  basis  (for  the  ideal  /  =  (G))  if 
the  G-normal  form  of  every  polynomial  is  unique.  We  now  increase  under- 
standing of  this  definition  by  giving  several  alternative  characterizations. 

We  first  note  two  general  results  on  partial  ordering  relations.  Let  X 
be  any  set  and  — >  be  a  binary  relation  on  X  and  — ►  be  its  reflexive 
transitive  closure.  Let  g,g',h  G  X.  We  call  h  a  common  successor  of  g  and 
g'  if  g  — *h  and  g'  — *h.  We  say  the  relation  — >  is  Church-Rosser  if  for 
all  elements  /,  5,  ^  G  X  we  have: 
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/  — ►  g  and  /  — >  y*  implies  that  g  and  g'  have  a  common  suc- 
cessor. 

The  relation  — ►  is  locall]/  Church-Rosstr  if  for  all  elements  f,g,g'  €  AT  we 
have: 

/  — ►  g  and  /  — ►  g*  implies  that  g  and  g'  have  a  conunon  suc- 
cessor. 

(Note:  "confluent"  is  an  alternative  terminology  for  "Church-Rosser".) 
The  relation  — >  is  Notthtrian  if  there  is  no  infinite  sequence  of  the  form 
/i  — »^  /j  — ►  •  •  •.  Note  that  if  — *  is  Noetherian  then  /  — *  f  must  not  hold 
for  any  /  G  AT.  An  element  /  6  X  such  that  /  — ►  g  implies  that  f  —  g 
is  said  to  be  in  normal  form.  Thus  normal  form  elements  are  "minimal" 
elements  (assuming  the  reduct  of  an  element  is  smaller  than  the  original). 
We  say  ff  is  a  normal  form  of  /  if  /  — ►  g\  the  set  of  all  normal  forms  of  / 
is  denoted  NF(/). 

Following  are  two  general  results  about  Noetherian  relations.  The  first 
proof  uses  a  very  powerful  principle  called  the  Principle  of  Noetherian  In- 
duction: 

For  a  Noetherian  relation  — y,  to  establish  the  validity  of  a 
predicate  P{x)  for  all  x  6  X,  it  is  sufficient  to  show:  if  P{y) 
holds  for  all  y  such  that  x  — >  y  and  x  j^  y  then  P{x)  holds. 

Let  us  demonstrate  the  validity  of  this  principle.  Suppose  that  it  is  shown 
that  P{x)  holds  whenever  P{y)  holds  for  all  y  7^  i  where  x  — »y.  For  the 
sake  of  contradiction,  suppose  P(io)  fails.  Then  for  some  Xi  ^  xo,  Xo  — ^  Xi 
and  P{xi)  fails.  Continuing  in  this  fashion,  we  obtain  an  infinite  descending 
sequence  xq,  xi,  Xj, . . . ,  which  contradicts  the  Noetherian  property. 

Lemma  10  Let  — >  be  a  Noetherian  relation.  Then  it  is  Church-Rosser  if 
and  only  if  it  is  locally  Church-Rosser. 

Proof.  One  direction  is  immediate.  In  the  other  direction,  assume  that 
the  relation  is  locally  Church-Rosser.  Let  /  — ►  g  and  /  — ►  h.  \{  f  —  g 
OT  f  =  h  then  the  result  is  irmnediate.    Otherwise,  let  /  — >  g'  — >  g  and 
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/  — >  h'  — ►  h.  Since  the  relation  is  locally  Church-Rosser,  g'  and  h'  have  a 
comnion  successor  /'.  By  the  principle  of  Noetherian  induction  (as  applied 
to  &')'  9  ^^^  f  have  a  common  successor  g";  similarly,  h  and  /'  have  a 
common  successor  h".  By  a  third  application  of  the  principle  to  /',  we 
get  that  g"  and  h"  have  a  common  successor  /".  Clearly  /"  is  a  common 
successor  of  §  and  h.  Q.E.D. 

Exercise.  Reprove  this  lemma  avoiding  the  Principle  of  Noetherian 
induction. 

Lemma  11  Let  — >  be  a  Noetherian  relation  on  X.  Then  — >  is  Church- 
Rosser  if  and  only  if  every  element  /   G  AT  has  a  unique  normal  form., 

|NF(/)|  =  1. 

Proof.  (=>)  By  our  assumption  that  the  relation  is  Noetherian,  NF(/)  is 
non-empty  for  any  /.  li  g,g'  G  NF(/),  then  the  Church-Rosser  property 
implies  that  g  and  g'  have  a  common  successor  h.  But  g  is  in  normal  fonn 
implies  g  =  h  and  similarly  g'  =  h. 

(<=)  Suppose  /-^ff  and  /-^ff'.  Then  NF(<7)  C  NF(/).  Since  NF(ff) 
is  non-empty  (Noetherian  assumption)  and  |NF(/)|  =  1,  we  conclude  that 
NF(g)  =  NF(/).  Thus  /  and  g  have  a  unique  common  normal  form,  say  h. 
Similarly  for  /  and  g'.  Thus  h  serves  as  the  common  successor  of  g  and  g'. 
Q.E.D. 

We  now  return  to  our  original  reduction  relation  — ►  (for  some  fixed 
finite  set  F  of  polynomials).  An  immediate  corollary  of  the  last  lemma 
gives  lis  the  first  equivalent  formulation  of  a  Grobner  basis. 

Corollary  1  G  is  a  Grobner  basis  if  and  only  if  the  relation  — ►  is  Church- 
Rosser. 

The  following  is  a  useful  tool. 

Lemma  12  If  f  —  g  — >  0  then  f  and  g  have  a  common  successor. 

Proof.  We  use  induction  on  the  number  of  steps  to  get  from  f  —  g  to  0. 
If  the  number  of  steps  is  zero  then  the  result  is  immediate.  Otherwise, 
suppose 

{f-g)-^h-^0. 
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It  18  suflBcient  to  show  that  h=  f  -ff  lor  some  /'  and  g'  such  that  /  — >  f 
and  g  -^  g*.  For  then,  the  induction  hypothesis  shows  that  /'  and  g'  have 
a  common  successor  which  is  also  the  common  successor  of  /  and  g.  Note 
that  h  =  [f  -  g)  -  atfo  for  some  fo  e  F,  a  ^  K  and  some  term  f  G  PP. 
Let  to  =  Hterm(</o)  be  the  eliminated  term.  It  is  clear  that  to  occurs  in  / 
and  g  with  some  (possibly  zero)  coefficients  Qi  and  02  (respectively)  such 
that  a  =  ai  —  ttj.  Thus 

h  =  {f-g)-{ai-  a^Yfo  =  {f-  aitfo)  -  {g  -  a^tfo) 

and  we  may  choose  /'  =  aitfo,  g'  =  oj^/o-  Q-E.D. 

Let  < >  be  the  relation  obtained  as  the  union  of >  and  its  reverse: 

/  <-^  g  if  and  only  if  either  /  — >  g  or  g  — >  /  holds. 

Let  *-^  denote  the  reflexive  transitive  closure  of  < — *. 
The  following  is  from  [Bachmair  and  Buchberger  19??] 

F 

Lemma  13  f  -  g  ^  [F)  if  and  only  if  f  *-^  g. 

Proof.  {<=)  This  is  easily  shown  by  induction  on  the  number  of  steps 
between  /  and  g.  Let 

/  =  ffo  " — ^  9\ ' — '  •  •  • " — *  9k  =  9 

for  some  k  >  0.  The  result  is  trivial  for  k  =  0.  Otherwise,  by  induction, 
9\-  9k^  [F)  and  it  is  seen  directly  from  the  definition  that  go  -  Qi  ^  [F). 
Thusgo-gte  {F). 

(=>)  If  f  -  g  e  {F)  then  we  can  express  f  -  g  a£  E.'^i  oaf,  where  each 
Q,  is  a  non-zero  monomial  and  the  fi  are  members  of  F.  The  /,'s  are 
not  necessarily  distinct  here.  If  m  =  0  the  result  is  trivial.  If  m  >  1 
then  we  can  write  g'  =  g  +  amfm  and  f  -  g'  =  HZV  "•/«•  By  induction 
hypothesis,  /  *^  g'.  We  also  have  that  g'  -  g  =  cimfm  — >  0.  By  the 
previous  lemma,  we  conclude  that  g  and  g'  have  a  common  successor  h. 
This  implies  g'  *-^  h  *-^  g.  This  shows  f  ^^  g'  ^  g.  Q.E.D. 
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Lemiaa  14  //  /  * — ►  g  then  f  and  g  have  a  common  successor. 

Proof.  We  use  induction  on  the  number  of  steps  to  get  from  /  to  g.  The 
result  is  trivial  for  0  or  1  step.    Suppose  /  •< — >  /'  •< — ►  g.    By  induction 

hypothesis,  /'  and  g  have  a  common  successor  h.    If  in  fact  / >  /'  then 

/  and  g  have  the  common  successor  h.  Otherwise,  suppose  /'  — ►  /.  But, 
since  /  and  h  are  both  successors  of  /',  an  application  of  the  Church-Rosser 
property  for  — >  implies  that  /  and  h  have  a  common  successor.    Then 

g  — *  h  implies  that  /  and  g  have  the  same  common  successor.  Q.E.D. 

Remark:  The  preceding  result  holds  for  any  relation  — y  that  is  Noethe- 
rian  and  Church-Rosser. 

Theorem  15  [Buchbergerj  G  is  a  Grobner  basis  if  and  only  if  for  all  f  G 
(G),  OeNFcl/). 

Proof.  {<=)  Let  h  be  any  polynomial.  We  want  to  show  that  |NFg(/i)|  =  1. 
Let  g,g'  G  NFg(/i).  Then  h  -  g  e  [G)  and  h  -  g'  e  {G).  Hence  g  -  g'  = 
[h  —  ^)  —  [h-  g)  €:  [G).  It  follows  that  g  —  g'  — *  0.  But  g  —  g'  \s  in  normal 
form  (since  g  and  g'  are).  It  follows  that  g  —  ^  =  0  or  g  =  g\^s  desired. 
(=>►)  Suppose  f  e{G).  Since  0  G  (G),  we  have  /-O  e  (G).  Then  lemma 

G 

13  implies  /•< — *0.  The  preceding  lemma  next  implies  that  /  and  0  have 
a  common  successor  h.  Since  0  is  in  normal  form,  /i  =  0.  This  proves  that 
OgNFg(/).  Q.E.D. 

Theorem  16  [Buchbergerj  G  is  a  Grobner  basis  if  and  only  if  for  all  f,g^ 
G,  0GNFG(5(/,ff)). 

Proof.  (=>•)  This  is  immediate  since  S{f,g)  G  (G)  if  f,g  G  G. 

{<=)   By  preceding  results,  it  is  enough  to  show  that  — ►  is  locally 

Church-Rosser.  So  let  / — yg  and  / — ► /i  and  we  are  required  to  show 
that  g  and  h  have  a  common  successor.  So 

g  =  f  -  oi/i  and  h  =  f  -  02/2 
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where  ai,ai  are  monomials  and  fi^fi  6  G.  Let  the  eliminated  monomials 
be  mi  =  Hmono(Qi/i)  and  mj  =  Hmono(a2/j).  There  are  two  cases: 

Case  mi  >m2.    (By  symmetry,  this  covers  the  case  mj  >  mi  as  well.) 

Then  we  see  that  mi  still  occurs  in  h.  Therefore  we  have  h  — >  h'  where 

h'  =  h-  aji  =  /  -  (qi/i  +  02/2). 

Now  observe  that  g  —  h'  =  0:2/2  — *^-  By  lemma  12  we  conclude  that  g 
and  h'  have  a  common  successor.  Then  g  and  h  have  the  same  common 
successor. 

Case  mi  =  mj.  Here  is  where  the  5-polynomials  come  in.  By  definition, 

where 

Hmono(/?i/i)  =  m  =  Hmono(^2/2)  (2) 

for  some  monomial  m  whose  corresponding  power  product  is  LCM(Hterm(/i),  Hlerm(/j 
Since  both  Hterm(/i)  and  Hterm(/2)  divide  mi  (=  m2),  we  conclude  that 
m  divides  mi.  Therefore  let  us  express  mi  as 

mi  =  mm 

for  some  monomial  m'.  Since  Hmono(ai/,)  =  mm'  and  Hmono(/i)  divides 
m,  we  have  that  m'  divides  a^.  Again,  we  may  factor  a^  as 

a,  =  m'^, 

for  some  monomial  7,  (t  =  1,2).  It  is  now  seen  that 

g  -  h    =    02/2  -  Qi/i 

=    rn'')2ft  -  mSi/i 

=     —[lih  -  lifi) 
m 

Since  H mono ( a, /j)  =  m,,  it  follows  from  the  last  equation  that  Hmono(2i^)  = 
1  or  m  =  Hmono('y,/,).  Comparing  this  with  (2),  we  conclude  that  -y,  =  /9,. 
In  otherwords, 

ff-;i  =  m'-5(/i,/2). 
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Since  we  assume  5(/i,/j)  — ►©,  lemma  12  implies  g  and  h  have  a  common 
successor.  This  concludes  our  proof  of  this  case. 

Note  that  the  above  two  cases  are  exhaustive:  in  particular,  it  is  not  pos- 
sible that  Hterm(mi)  =  Hterm(m2)  and  Hcoef(mi)  ^  Hcoef(m2).  Q.E.D. 

Remark:  Compared  to  the  definition  of  a  Grobner  basis,  theorenis  15 
and  16  are  successively  better  characterizations  of  a  Grobner  ba^is.  More 
precisely,  in  the  original  definition  we  insist  that  every  polynomial  must 
have  a  unique  normal  form,  theorem  15  has  a  similar  requirement  but  only 
for  polynomials  in  the  ideal,  and  theorem  16  has  the  same  test  restricted  to 
only  a  finite  set  of  polynomials.  Thus,  theorem  16  is  an  effective  criterion  to 
test  if  a  given  F  is  indeed  a  Grobner  basis:  given  F,  for  each  pair  /,  y  G  F, 
we  compute  idf{S{f,g))  and  see  if  it  is  equal  to  zero.  If  this  is  zero  for  all 
/,  g  then  we  know  F  is  a  Grobner  basis.  Conversely,  if  any  nf/?(5(/,  g))  ^t  0 
then  either  |NFj.(5(/,s))|  >  0  or  0  ^  NFj.(5{/,y)).  In  either  case  F  is  not 
a  Grobner  basis.  Of  course,  this  still  does  not  tell  us  how  to  construct  a 
Grobner  basis  for  the  ideal  generated  by  a  given  F.  The  next  section  will 
do  this. 

To  conclude  this  section,  we  now  give  characterizations  of  Grobner  basis 
ba^ed  on  looking  at  headterms. 

Theorem  17  G  =  {</i, . . . ,  gr}  is  a  Grobner  basis  if  and  only  if  every  poly- 
nomial f  ^  [G)  can  be  expressed  in  the  form 

T 

f  =  Y1  <*•&• 

»=i 

where  Hterm(a,<7,)  <  Hterm(/)  for  all  i. 
Proof.  {=>)  Let  /  e  {G)  and  let 

be  a  reduction  sequence  witnessing  the  fact  that  /  — ►  0.  So  for  j  =  1, . . . ,  it, 
hj  6  G  and  fj  =  /y_i  —  mjhj  for  some  monomial  m^.    We  say  that  rrij 
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is  associated  with  hj.  Then  we  see  that  /  =  J^[_j  q,<7,  where  p,  6  G 
and  the  a,  are  obtained  by  summing  all  the  monomials  associated  with 
y,.  Now  Hterm(my/i;)  <  Hterm(/y_i)  and  Hterm(/,)  <  Hterm(/j_i)  for  all  j 

and  hence  Hterm(my/ij)  <  Hterm(/o).  We  conclude  that  Hterm(a,g,)  <  Hterm(/o), 
as  we  wanted  to  show. 

(<=)  Let  /o  €  (G).  By  assumption,  /o  =  X)i=i  ct,J?i  where  Hterm(Q,g,)  <  Hterm(/o) 

A 

This  implies  that  there  is  some  set  of  indices  J  C  {l,...,r}  such  that 
Hmono(/o)  =  X^yeJ  Hmono(aj(7y).  Choose  any  j  €  J  and  it  is  not  hard  to 
see  that  the  head  monomial  of  /o  can  be  eliminated  by  applying  Qj.  Let  /i 
be  the  result  of  this  elimination.  Since  /i  G  (G),  we  may  apply  the  same 
argument  to  eliminate  the  head  monomial  of  /i  to  obtain  /2.    Repeating 

this,  we  get  the  sequence  /o ► /i — > /z — ►•••.    This  sequence  must  be 

finite  and  it  is  easy  to  see  that  the  last  term  is  0.  This  shows  that  G  is  a 
Grobner  basis.  Q.E.D. 

Theorem  18  Let  G  =  {gi,...  ,gr}-  Every  polynomial  f  G  (G)  can  be 
expressed  in  the  form 

r 

/  =  Z]  OCiQi 

»=1 

where  Hterm(/)  >  Hterm(Q,j?,)  for  alii  if  and  only  if  the  set  G'  =  {Hmono(y)  : 
g  G  G}  16  a  basis  of  the  ideal  generated  by  P  =  {Hmono(/)   :  /  G  (G)}. 

Proof.  (=>•)  Suppose  /  G  (/'),  we  have  to  express  it  as  2m  element  of  (G'). 
We  can  write  /  in  the  form  /  =  Z);=i  Hmono(/,)  for  some  k  and  fj  G 
(G).  But  each  fj  can  be  expressed  as  /y  =  Z)--i  Q:;,,ffi  where  Hterm(/y)  > 
Hterm(Qj  ,5,  for  all  t.  Then  we  see  that 

k  T 

f  ^J2  Hmono(  J3  <^j.i9.) 

This  implies  that  /  can  be  generated  from  G'  =  {Hmono(g,)  :  t  =  1, . . . ,  r}. 

(«^)  Given  /  G  (G),  by  assumption,  Hmono(/)  can  be  expressed  as 

Z)i=i  »n,Hmono(<7,)  for  suitable  monomials  mj.     Li  fact  we  see  that  all 

but  one  of  the  m.'s  may  be  assumed  to  be  zero.    Let  /o  =  /  and  /i  = 
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/o  -  Ei=i  rriig,.  Obeerve  that  Htenn(/o)  >  Hterm(/i).  Since  /i  e  (G),  we 
can  repeat  this  argument  to  obtain  /j  with  Hterm(/i)  >  Hterm(/2).  This 
process  must  halt  after  k  repetitions;  the  last  polynomial  ft  is  seen  to  be 
0.  It  follows  that  /o  can  be  expressed  as  a  sura  of  the  form  12^=1  <*!?•  where 
Hterm(Qif?,)  <  Hterm(/o).  Q.E.D. 

A 

Remark:  The  characterization  of  Grobner  basis  in  this  last  theorem  is 
used  in  the  literature  to  define  the  concept  of  a  standard  basis. 

6      The  Basic  Algorithm  of  Buchberger 

Now  we  are  ready  to  present  the  basic  form  of  the  Grobner  ba^is  algorithm 
of  Buchberger.  The  algorithm  is  deceptively  simple  -  its  complexity  analysis 
is  not.  Like  the  nornaal  form  algorithm,  it  is  non-deterministic.  In  fact,  it 
calls  the  normal  form  algorithm  at  a  crucial  point. 

Input:  F  a  finite  set  of  polynomials. 
Output:  G  a  Grobner  basis  for  F. 

G:=F 

B:={{f,9}  :  LgeGJ^g} 

while  B  7^  0  do  begin 

Choose  {f,g}  to  be  any  pair  in  B 

B:=B-{{f,g}} 

h:=S{f,g) 

h':=-nlG{h) 

if  h'  ^  0  then  begin 
B:^Bu{{f',h'}  :  f'eG} 
G:=Gu{h'} 

end  {if} 
end  {while} 


Lemma  19   The  basic  algorithm  terminates. 
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Proof.  Since  we  have  already  shown  that  the  normal  form  algorithm  termi- 
nates, it  suffices  to  show  that  the  number  of  iterations  of  the  while-loop  is 
finite.  Let  hi  be  the  polynomial  assigned  to  the  program  variable  h'  in  the 
ith  iteration  of  the  loop.  For  the  sake  of  contradiction,  assume  that  there 
are  infinitely  many  iterations.  It  is  seen  that  in  each  iteration  in  which  h' 
is  zero  the  size  of  the  set  B  decreases.  Since  B  is  always  finite  in  size,  there 
cannot  be  an  infinite  succession  of  iterations  in  which  h'  =  0.  Hence  there 
are  infinitely  many  non-zero  values  of  /i,. 

Now  let  H  be  the  infinite  set  consisting  of  all  the  non-zero  /i,'s.  By 
Dixon's  lemma,  there  is  a  finite  set  of  indices  /  such  that  each  g  ^  H  is  a. 
multiple  of  /i,  for  some  i  G  /.  Let  t'o  be  the  largest  index  in  /.  Then  for 
all  j  >  t'o,  if  hj  ^  0  then  hj  is  a  multiple  of  hi  for  some  i  G  /.  However, 
during  the  jih  interation,  the  polynomial  h,  is  in  the  set  G.  Hence  hj  is 
not  a  G-normal  form,  contradiction.  Q.E.D. 

Theorem  20   The  Basic  Algorithm  is  correct. 

Proof.  Since  the  algorithm  terminates,  we  only  have  to  show  that  at  termi- 
nation the  set  G  is  a  Grobner  basis.  By  a  characterization  in  the  previous 
section,  we  must  show  that  for  all  /,  <?  G  G,  S{f,g)  — ►O.  The  loop  in- 
variant that  we  may  observe  is  this:  for  all  f,g  e  G,  if  {f,g}  ^  B  then 
S{f,g)  -^0.  This  invariant  holds  at  the  beginning  of  the  iterations.  Each 
iteration  preserves  the  invariant.  Since  the  set  B  is  empty  on  termination, 
the  set  G  is  indeed  a  Grobner  basis.  Q.E.D. 

Unfortunately,  there  are  no  upper  bounds  on  the  rimning  time  of  the 
algorithm  except  for  some  special  c«ises.  Unless  otherwise  noted,  complex- 
ity of  this  algorithm  will  be  just  the  worst  case  number  of  iterations  of  the 
while-loop.  For  instance,  upper  bounds  are  known  for  the  two  and  three 
variable  cases  and  assuming  a  total  degree  ordering  of  terms;  the  bounds  in 
these  cases  are  a  polynomial  and  a  single  exponential,  respectively.  Other 
bounds  are  known  in  terms  of  more  technical  quantities  such  as  the  'reg- 
ularity' of  ideal.  As  for  lower  bounds,  there  is  a  double  exponential  lower 
bound  essentially  due  to  Mayr  and  Meyer. 
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7      Uniqueness  of  Reduced  Grobner  Bases 

In  this  section  we  show  that  a  Grobner  basis  that  is  "reduced"  in  a  certain 
sense  is  unique  for  an  ideal.  We  give  a  method  to  obtain  such  bases. 

We  call  a  basis  F  minimal  if  for  all  f  e  F,  {F)  /  (F  -  {/}).  A  basis  F 
is  self-reduced  if  for  all  /  G  F,  /  is  a  (F  -  {/})-nonnal  form. 

We  call  a  Grobner  basis  G  reduced  if 

1.  for  all  <7  e  G,  Hcoef(<7)  =  1 

2.  G  is  minimal 

3.  G  is  self-reduced 

Note  that  if  F  is  a  minimal  bcisis  and  (F)  ^  (0)  then  0  ^  F. 

Lemma  21  //  g,g'  G  G  are  distinct  polynomials  in  a  minimal  Grobner 
basis  G  then  Hterm(y)  does  not  divide  Hterm(<7').  In  particular,  no  two 
polynomials  in  G  have  the  same  headterm. 

Proof.  Note  that  (G)  7^  (0)  so  g,  g'  are  distinct  from  0.  Let  Hterm(y)  divide 
Hterm(^')  and  G'  =  G-{g'}.  Since  G  is  minimal,  (G')  7^  (G).  So  g'  ^  (G') 
and  if  /i  e  NFc'ig')  then  h  ^  0.  But  since  G  is  a  Grobner  basis,  h  — ►  0.  In 
particular,  k  is  reducible  by  some  polynomial  f  €:  G.  If  /  =  </'  then  since 
Hterm(ff)  divides  B.ienn{g'),  h  is  also  reducible  by  g.  Since  g  G  G',  this 
contradicts  the  assimaption  that  /i  is  a  G'-normal  form.  On  the  other  hand, 
\i  f  ^  g'  then  f  G  G  and  we  arrive  at  the  same  contradiction.  Q.E.D. 

Lemma  22  //  G,  G'  are  minimal  Grobner  basis  for  the  same  ideal,  (G)  = 
(G'),  then  the  set  of  headterms  in  G  is  equal  to  the  set  of  headterms  in  G'. 

Proof.  Suppose  for  some  g  ^  G,  Hterm(ff)  does  not  occur  aniong  the  head- 
terms  of  polynomials  in  G'.  Since  g  — ►  0  and  g  ^  0,  there  is  some  g'  G  G' 

such  that  Hterm(g')  properly  divides  Hterm((7).  Again,  since  g'  — ►O,  there 
is  some  g"  G  G  such  that  Hterm(5'')  divides  (not  necessarily  properly) 
Hterm(g).  This  means  that  B.ienn{g'')  properly  divides  Hterm(y),  contra- 
dicting the  preceding  lemma.  Q.£.D. 
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Corollary  2  //  G,  G'  are  two  minimal  Grobntr  bases  for  the  same  ideal 
then  \G\  =  \G'\. 

Of  course  the  size  of  the  minimal  Grobner  bases  for  any  ideal  may  still 
depend  on  the  choice  of  admissible  orderings. 

Theorem  23  The  reduced  Grobner  basis  of  an  ideal  is  unique  (relative  to 
the  choice  of  an  admissible  ordering). 

Proof.  Let  G,G'  be  two  reduced  minimal  Grobner  bases  for  the  same  ideal. 
We  obtain  a  contradiction  by  supposing  that  there  is  some  polynomial  g  in 
G  —  G'.  By  the  previous  lemma,  there  is  some  other  polynomial  g'  in  G'  —  G 
such  that  Hmono(ff)  =  Hmono(g')  (recall  that  Hcoef(p)  =  1  =  Hcoef(ff')). 
Let  h  =  g  —  g'.  Then  h  ^  0  and  h  — >  0  since  G  is  a  Grobner  basis.  So  some 
term  t  occuring  in  h  can  be  eliminated  by  application  of  some  f  ^  G.  Now  t 
must  occur  in  g  or  g'.  If  t  occurs  in  g  then  g  is  reducible  by  /,  contradicting 
the  assumption  that  G  is  reduced.  If  t  occurs  in  g"  then  let  /'  G  G'  such 
that  Hterm(/')  =  Hterm(/).  Again  g'  is  reducible  by  /',  contradicting  the 
assumption  that  G"  is  reduced.  Q.E.D. 

First  we  present  a  method  to  compute  the  reduced  Grobner  basis  cor- 
responding to  any  input  set  F  of  polynomials: 
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Function  SelfReduce(F) 

Input:  F  a  finite  set  of  polynomials. 

Output:  R  a  reduced  beisis  for  (F). 

i2:=0 

while  F  7>^  0  do  begin 

Choose  /  from  F  and  set  F  :=  F  —  {/} 

h:=S{f,g) 

if  /  /  0  then  begin 

TEMP  :={{ge  R  :  g  is  reducible  by  /} 
iE  :=  (i?  -  TEMP)  U  {/} 
F:=FUTEMP 
end  {if} 
end  {while} 
return(i2) 


Lemma  24   The  routine  SelfReduee  terminates. 

Proof.  We  will  trace  the  history  of  transformations  for  a  single  polynomial 
/  from  the  original  input  set  F.  In  each  iteration  of  the  while-loop,  there 
is  a  selected  polynomial  /,  and  in  a  natural  way,  this  polynomial  can  be 
associated  with  a  polynomiaJ  in  the  original  input  set  F.  Observe  that 
except  possibly  for  the  first  time  that  the  polynomial  /  is  selected,  all  sub- 
sequent selection  of  (versions  of)  /  w^ll  result  in  a  further  transformation, 
i.e.,  in  the  loop,  /  ^  nialf).  If  /  =  /o,  /i,  /j,  •  •  • ,  constitute  the  successive 
transformed  versions  of  /  then  the  same  argument  iised  in  the  termination 
proof  for  the  Normal  Form  algorithm  implies  that  this  sequence  of  /,'s  is 
finite.  Since  this  is  true  for  every  /  in  the  original  F,  there  can  only  be  a 
finite  number  of  iterations.  Q.E.D. 

Lemma  25    The  routine  SelfReduee  is  correct. 
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Proof.  Two  loop-invariants  can  be  easily  checked  to  hold:  (i)  the  set  R  is 
self-reduced  and  (ii)  the  ideal  {F  U  R)  does  not  vary.  The  final  result  then 
has  the  desired  property.  Q.E.D. 

We  are  finally  ready  to  present  the  algorithm  for  computing  reduced 
Grobner  bases.  This  is  done  by  simple  modifications  to  the  Basic  Algo- 
rithm. 


Input:  F  a  finite  set  of  polynomials. 
Output:  G  a  Grobner  basis  for  F. 

G:=  SelfReduce(F) 
B:^{{f,g}  :  f,9eG,f^g} 
while  jB  /  0  do  begin 

Choose  {f,g}  to  be  any  pair  in  B 

B~B-{{f,g}} 

h:=S{f,g) 

h':=niG{h) 

if  /i'  7^  0  then  begin 
B:=BU  {{/'>'}  :  fee} 
G  :=  SelfReduce(G'  U  {h'}) 

end  {if} 
end  {while} 


We  leave  the  correctness  proof  for  this  procedure  to  the  reader.  Note 
that  the  application  of  SelfReduce(G  U  {h'})  within  the  while-loop  can  be 
made  more  efficient  if  we  exploit  the  fact  that  G  is  already  self-reduced. 
We  also  leave  this  improvement  of  the  procedure  SelfReduce  to  the  reader. 

Exercise.  Let  G  be  a  Grobner  basis  such  that  for  some  pair  g,g'  of 
distinct  polynomials  in  G,  Hterm{ff)  divides  Hterm(ff').  Then  G'  =  G  -  {g'} 
is  still  a  Grobner  basis.  Remark:  This  result  does  not  say  that  (G')  =  (G), 
and  indeed  this  is  false  in  general. 
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8      Applications 

In  this  section,  we  discuss  applications  of  the  Grobner  basis  algorithm  to 
decide  several  fundamental  problems  in  Ideal  Theory.  The  wide  range  of 
applicability  of  the  Grobner  basis  algorithm  is  a  clear  testimony  to  the 
power  of  this  concept. 

8.1      Ideal  Theoretic  Problems 

Let  R  =  K[xi,  . . .,  i„]  be  a  polynomial  ring  over  a  field  K,  and  (F)  C  R 
be  an  ideal  defined  by  a  finite  set  of  generators  F.  In  what  follows,  G  will 
always  stand  for  the  Grobner  baisis  of  [F). 

We  define,  in  the  usual  way,  an  equivalence  relation  for  any  ideal  I  C  R, 
=  mod  /,  {congruence  modulo  the  ideal  I)   over  R  a.s  follows:    for  all  /, 

f  =  g  mod  I  iS  f  -  g  e  I. 

Let  F  C  R,  and  G  the  Grobner  basis  of  (F).  It  is  easy  to  see  that  if 
/  — >  g  then  f  =  g  mod  (G).  By  transitivity,  if  /  — >  g  then  f  =  g  mod  (G) 
In  particular,  /  =  nfcf/)  mod  (G). 

Next,  if  nfG(/)  =  nfc{g)  then,  by  the  above  argument,  we  conclude  that 
f  =  g  mod  (G).  Conversely,  if  f  =  g  mod  (F)  then,  by  Lenmaas  13  and  14, 
we  see  that  /  and  g  have  a  common  successor  h  with  respect  — k  Since  G 
is  Grobner,  a  corollary  of  Lemma  11  implies  that  nfG(/)  =  rdc{h)  =  nfG(y). 
Thus 

/  =  ff  mod  (F)  ^  nfcif)  =  nicig). 

Let  T  be  an  arbitrary  set  with  an  equivalence  relation  ~  defined  on 
it.  For  instance,  for  our  application,  we  may  choose  T  to  he  R  —  K[xi, 
...,  Xn]  and  ~  to  be  the  equivalence  relation,  =  mod(F).  We  will  also 
assume  that  T"  is  a  decidable  set,  i.e.  there  exists  a  decision  procedure  for 
the  membership  problem  for  T.  This  will  be  true  of  all  the  polynomial  rings 
of  interest  to  us. 

Following  Buchberger,  we  define  a  canonical  simplifier  for  ~  on  T  to  be 
an  algorithm  C  with  input  and  output  in  T  such  that  for  all  f,gGT 

•  /  ~  C{f)  and 
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•  f-g^C{f)  =  C{g). 

Notice  that  the  function  C  gives  a  unique  representative  in  each  equivalence 
class  of  T/  ~.  C{f)  is  called  a  canonical  form  of  /. 

Problem  1   /CANONICAL  SIMPLIFIER/ 

•  INPUT:   F  C  R. 

•  OUTPUT:  A  canonical  simplifier  algorithm  C  for  =   raod  [F)  on  K[xi, 

First,  compute  G  the  Grobner  basis  of  F.  For  any  /  let  C{f)  =  nioif)- 
The  correctness  of  this  method  follows  directly  from  the  discussion  above. 

Problem  2  /IDEAL  CONGRUENCE/ 

•  INPUT:  F  CR,  and  f,  geR. 

•  DECIDE:  f  =  gmod{F)? 

As  before,  compute  G  the  Grobner  basis  of  {F).  Determine  n(o{f)  and 
iiio{g).  Output  true  if  nf<;(/)  =  nf<;(g);  otherwise,  false.  The  correctness 
of  this  method  is  follows  from  the  preliminary  observations  made  above. 

Problem  3  /IDEAL  MEMBERSHIP/ 

•  INPUT:   F  CR,  and  f  e  R. 

•  DECIDE:   f  e  [F)? 

We  note  that  van  der  Waerden  calls  this  'the  central  problem  of  ideal 
theory.' 

Compute  G  the  Grobner  basis  of  F,  and  the  normal  form  of  /,  nfG(/), 
with  respect  to  -^.  Output  true  if  nf(?(/)  =  0;  otherwise,  false.  To  see 
the  correctness  of  the  algorithm, 

fe{F)   iff  fe{G) 

iff    nfc(/)  =  0     (from  Theorem  15) 
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Problem  4  /SUBIDEALS/ 

•  INPUT:  Fi,  FjC  R. 

•  DECIDE:   (Fi)  C  {F,)  (Is  (Fj)  a  suhidtal  of  [F2))? 

Compute  Gj  the  Grobner  basis  of  {Fi).  Output  true  if  for  all  f  ^  Fi, 
nfcjf/)  =  0;  otherwise,  false.  It  is  easily  seen  that 

{Fi)C{F2)<^^feFje{F2). 

Problem  5   /IDEAL   EQUALITY/ 

•  INPUT:  Fi,  F2C  R. 

•  DECIDE:    (Fi)  =  (F2)? 

This  problem  clearly  reduces  to  the  previous  problem. 

Alternatively,  compute  G[  and  Gj  the  reduced  Grobner  bases  for  (Fj) 
and  (Fj),  respectively.  Output  true  if  G[  —  G\\  otherwise,  false.  By 
Theorem  23,  (Fi)  =  (Fj)  if  and  only  if  G\  =  G,,  and  the  correctness 
follows. 

8.2      Residue  Class  Ring  Modulo  an  Ideal 

Let  R  =  K\xx,  . . .,  x„]  be  a  polynomial  ring  over  the  field  K  and  /  an 
ideal  of  R.  Then  the  equivalence  relation  =  mod  /  partitions  the  ring  R 
into  equivalence  classes  such  that  f ,  g  ^  R  belongs  to  the  same  class,  if 
f  =  g  mod  /.  The  equivalaince  classes  of  R  are  called  its  residue  classes 
modulo  /.  We  use  the  notation  R/I  to  represent  the  set  of  residue  classes 
of  R  with  respect  to  /.  Let  /  denote  the  set  {g  \  f  =  g  mod  /}.  It  is  easy 
to  see  that  the  map  /  •— ►  /  is  the  natural  ring  homomorphism  of  R  onto 
R/I.  We  sometimes  write  /  +  /  for  /.  R]I  is  called  the  residue  class  ring 
modulo  I. 

Lemma  26  Rj I  is  a  vector  space  over  K. 
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Lemma  27  Let  G  be  a  Grobner  basis,  and  let 

B    =    {p  I  P  €  PP(xi, . . . ,  J„)  such  that  p  is  not  a  multiple 

of  the  Htenn  of  any  of  the  polynomials  in  G}.  (S) 

Then  B  is  a  linearly  independent  (vector  space)  basis  of  Rj{G)  over  K. 

Proof.  Let  7  =  /  +  (G)  be  an  element  of  Rl[G)  and  let 

nfcC/)    =    ci -Pi  +  •••  +  <;« -Pi, 

where  Ci  e  K  and  pi  6  PP(ii,..  .,x„). 

Hence,  for  all  1  <  »  <  /,  p,  is  not  a  multiple  of  the  Hterm  of  any  polynomial 
in  G,  and  p,  G  B.  Since  we  can  write  /  as 

J  =  ci-  pi  + \-ci-  pi, 

B  spans  the  vector  space  R/{G). 

Furthermore,  the  elements  of  B  are  linearly  independent.  Assume  to 
the  contrary,  that  is,  for  some  pi,  . . .,  Pm  6  B,  we  can  obtain  Ci,  . . .,  Cm, 
not  all  0,  such  that 

Ci'  pi  + 1- Cm- Pin  =  0. 

In  other  words, 

/  =  Cl-pi  +  -"  +  Cm-Pm€(G). 

Hence  nfcif)  =  /  =  0,  i.e.,  Ci  =  •  ••  =  Cm  =  0,  thus  resulting  in  a  contra- 
diction. Q.E.D. 

Problem  6  /FINITE   DIMENSIONALITY  OF  A   RESIDUE   CLASS   RING/ 

•  INPUT:   F  C  R. 

•  DECIDE:   Is  R/{F)  finite-dimensional? 

Compute  G  the  Grobner  basis  of  (F).  Output  true  if  for  all  t  (1  <  t  < 
n),  a  power  product  of  the  form  i-'  (j,  >  0)  occurs  among  the  Hterm's  of 
the  polynomials  in  G;  otherwise,  false.  To  see  this,  if  for  some  t,  none  of 
the  power  products  of  Hterms  in  G  has  the  required  form  then  we  get  an 
infinity  of  basis  elements.  The  converse  is  similar. 
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Problem  7  /BASIS  OF  A  RESIDUE  CLASS  RING/ 

•  INPUT:  F  C  R. 

•  OUTPUT:  If  R/{F)  is  a  finite  dimensional  vector  space  then  output 

1.  A  basis  B  of  the  vector  space  R/{F). 

e.  The  'multiplication  table'  for  R/{F);  the  {p„Pj)th  entry  of  the 
table  gives  a  linear  representation  of  pi  •  pj  in  terms  of  the  basis 
elements  in  B. 

Compute  G  the  Grobner  basis  of  [F).     Let  the  set  B  be  as  in  the 
Equation  3.  This  is  easily  computed. 

For  each  Pi,  pj  ^  B,  compute  the  normal  form  of  p,  ■  Pj' 

nfclp.  •  P:)  =  ci  •  pi  +  •  •  •  +  c„  •  p„. 

Then  (^, Pj)th  entry  of  the  multiplication  table  is 

Ci-Pi  +  ■••  +  Crn-Pm- 

The  correctness  of  the  procedure  follows  immediately  from  the  Lenoma  27. 

8.3      Solving  Systems  of  Polynomial  Equations 

Grobner  bases  can  be  advantageously  used  to  settle  several  important  ques- 
tions about  solvability,  number  of  zeroes  and  finally,  finding  the  zeroes  of  a 
system  of  polynomials  F.  We  first  show  a  result  that  is  essentially  equiva- 
lent to  Hubert's  Nullstellensatz. 

Lemma  28  Every  maximal  ideal  M ,  distinct  from  the  unit  ideal  R  has  a 
zero  in  an  algebraic  extension  of  K. 

Proof.  To  every  polynomial  /(ii,  X2>  •  •  •»  ^n)  assign  an  element  of  the 
residue  class  ring  R/M,  given  by  the  natural  ring  homomorphism.  Since 
M  ^  R,  every  element  a  of  K  will  correspond  to  the  distinct  element 
a  =  a  +  M.  (Otherwise,  li  a  ^  b  and  a  =  b  then  a  —  6eMsol  = 
(a  —  6) (a  —  6)~*  6  Af  and  then  M  =  R.)  Since  M  is  a  maximal  ideal,  R/M 
is  an  algebraic  (field)  extension  of  K. 
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Let  xi,  Xi,  . . .,  Xn  be  the  images  of  Zi,  Zj,  . . .,  z„  under  the  natural  ring 
homomorphism  from  R  into  R/M. 

Since,  the  ring  operations  in  R/M  are  naturally  induced  from  the  same 
operations  in  R,  and  ae  K  maps  into  itself,  we  see  that,  for  every  f  e  M, 
f{xi,  Zj,  ...,  Xn)  =  0  in  R/M.  Hence,  M  has  a  zero  in  an  extension  field 
of  X.  Q.E.D. 

Theorem  29  F  is  solvable  if  and  only  if  1  ^  (F). 

Proof.  If  1  e  [F)  then  it  is  easily  seen  that  F  is  unsolvable.  We  prove 
the  converse  by  contradiction.  Among  all  the  ideals  I  such  that  1^7,  but 
I  is  unsolvable,  there  is  a  maximal  ideal  M  (follows  from  the  maximum 
principle).  Since  1  ^  M,  M  is  distinct  from  the  unit  ideal.  Since  M  is 
maximal,  from  the  previous  Lemma  (28),  we  conclude  that  M  has  a  zero 
in  some  extension  field  of  K,  thus  deriving  a  contradiction.  Q.E.D. 

Problem  8  /SOLVABILITY/ 

•  INPUT:  F  C  R. 

•  DECIDE:  Does  F  have  a  solution  in  an  algebraic  extension  of  K? 

Compute  G  the  reduced  Grobner  basis  of  (F).  Output  true  if  1  ^  G; 
otherwise,  false.  From  Theorem  29,  it  follows  that  F  is  unsolvable  if  and 
only  if  1  €  (G).  But  1  G  (G)  if  and  only  if  nfG(l)  =  0,  that  is,  if  and  only 
if  1  e  G. 

We  state  the  following  useful  theorem  without  proof: 

Theorem  30  F  has  finitely  many  zeroes  (in  every  algebraic  extension  of 
K)  if  and  only  if  the  vector  space  R/{F)  has  a  finite  vector  space  dimension. 

Problem  9  /NUMBER  OF  ZEROES/ 

•  INPUT:   F  C  R. 

•  DECIDE:   Does  F  have  a  finitely  many  solutions? 
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Output  true  if  R/{F)  is  a  finite-dimensional  vector  space  (this  can  be 
found  out  using  the  solution  to  Problem  6);  otherwise  false.  The  correctness 
of  this  method  follows  from  the  Theorem  30. 

Next  we  discuss  the  problem  of  finding  all  the  real  solutions  of  a  systems 
of  polynomials  F.  We  assume  an  'oracle'  that  provides  all  the  zeroes  of  a 
univariate  polynomial  p.  In  general  we  cannot  represent  the  roots  using 
finite  precision  numbers.  However,  there  are  techniques  to  represent  alge- 
braic numbers  so  that  we  can  answer  any  question  about  the  roots  without 
error.  This  is  not  really  a  deficiency  of  Buchberger's  algorithm  but  an 
intrinsic  limitation  of  algorithmic  solvability  of  polynomial  equations. 

For  the  method  below,  we  assume  that  the  Grobner  bcisis  is  computed 
using  a  lexicographic  admissible  ordering  on  the  power  products;  the  method 
can  be  modified  to  deal  with  other  admissible  orderings.  For  simplicity,  we 
omit  discussions  of  the  other  methods.  The  method  discussed  below  is 
adapted  from  Buchberger's  survey  paper. 

The  following  useful  lemma  shows  that  the  '»''"  eliminatton  idea/'  of  a 
Grobner  basis  G  is  generated  by  just  those  polynomials  in  G  that  depends 
on  the  variables  ii,  . . .,  x,. 

Lemma  31  Let  G  be  Grobner  basis  with  respect  to  the  purely  lexicographic 
ordering  of  power  products.  Without  loss  of  generality  let  us  cissume  that 
^1  <i«x  ^2  <i«x  •  •  •  <i«  Xn-   Then  for  1  <  i  <  n 

{G)nK\xu...,Xi]  =  {GnK{x„...,Xi]). 
Proof  Let  /  €  (G)  n  K[xi, . . . ,  i,].  Then  /  -7^  0,  i.e. 

f  -  f  JL  f  A^       "^  f    A^  f        -  n 

J   —  JO ^  /l  >  •  •  •  ^  /m *  Jm+l  —  U, 

where  ^^'s  are  in  G.  We  claim  that,  for  all  J  (0  <  j  <  m  -h  1),  /  can  be 
written  as 

/  =  /,  -f- ai  •  ^1 -I- \-aigu 

where  /y  €  [G)  n  K\xi,  . . .,  x,|,  y^'s  are  in  G  fi  K\xi,  . . .,  x,]  and  a^'s  are 
in  K[xi,  ...,  X.]. 

Clearly  f  =  fo  satisfies  the  claim.   Consider  the  j^^  {j  >  0)  reduction: 

fj  —^  />+!•  Since,  by  the  inductive  hypothesis,  fj  G  {G)nK[xi,  . . .,  x,],  the 
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Hmono  of  g'^  mxist  involve  only  the  first  t  variables.  Furthermore,  since  the 
admissible  ordering  chosen  is  lexicographic,  every  monomiaJ  of  g'^  involves 
the  same  set  of  variables.  Hence,  fj+i  =  f,  —  m,  •  <j^,  and  rrij  is  a  monomial 
in  the  first  t  variables,  thus  the  claim  follows  immediately. 
In  particular,  we  can  write  /  as 

/  =  ai-<7iH \-(Xi-gu 

where  g^s  are  in  G  n  K\xi,  . . .,  x,]  and  q^'s  are  in  K\xi,  ...,  ij].   Hence 
f  ^  {G  n  K\xi,  . . .,  I,]),  where  the  ideal  is  formed  in  K[xi,  . . .,  x,]. 
This  shows  that 

(G)  n  K[ii, . . . ,x,]  c  (G  n  K[xi,. . . , X.]). 

The  converse  is  trivial: 

{G)nK{xu...,Xi]D{GnK[xu...,Xi]). 

Q.E.D. 

Problem  10  /FINDING  ALL  THE  ZEROES  OF  F] 

•  INPUT:  F  C  R.  F  is  solvable,  with  finitely  many  solutions. 

•  OUTPUT:  All  the  solutions  of  the  system  F. 

Compute  G  the  Grobner  basis  of  (F)  with  respect  to  the  purely  lexi- 
cographic ordering  of  the  monomials.  Let  G,  (1  <  i  <  n)  be  defined  as 
follows: 

Gi    =    Gn^lxi] 

G.     =     GnK[xi,...,x.]-ii:[xi,...,x,_i],      2<j<n 

By  the  previous  lemma  (G)  n  K[xi,  . . .,  x,]  =  (Uy^jGy).    Furthermore  G 
contains  exactly  one  polynomial  p  of  i^[xi]. 

The  successive  elimination  can  be  carried  out  by  the  following  algorithm: 
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p:=  the  polynomial  in  Gi 

Xr  :=  {(a)  |  p{a)  =  0} 

for  t  :=  1  to  n  —  1  do  begin 

forall  (oi, . . .  ,a,)  €  AT,  do  begin 

H:={g{ai,...,ai,Xi+i)  \  ge  G.+i} 

p:=  the  greatest  common  divisor  of  the  polynomials  in  H 

if  p  7^  1  then  begin 

Ar.+  i:=X.+  iU  {(ai,...,a.,a)  |  p(a)  =0} 
end  {if} 
end  {forall} 
end  {for} 
return  X„ 
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